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A continuous train of irregularly spaced spikes, peculiar of homoclinic 
chaos, transforms into clusters of regularly spaced spikes, with quiescent pe- 
riods in between (bursting regime), by feeding back a low frequency portion 
of the dynamical output. Such autonomous bursting results to be extremely 
robust against noise; we provide experimental evidence of it in a CO2 laser 
with feedback. The phenomen here presented display qualitative analogies 
with bursting phenomena in neurons. 

PACS numbers : 05.45.-a; 89.20.-a; 89.75.-k. 



Spiking occurs in many physical and biological systems in neurons it is associated with 
the generation of action potentials [|J. In general, spikes occur irregularly and the precise 
time position of each spike within a train is considered to code relevant information Q. A 
single neuron can also show a more complicated firing pattern, that is, bursting. The term 
bursting refers to short trains of rapid spike oscillations intercalated by quiescent intervals, 
which periodically repeat. This behavior is found in many excitable biological systems as well 
as in chemical reactions |lp|J • In the neural context bursting phenomena have been found in 
different cell types: examples include thalamic neurons J/J, hippocampal pyramidal neurons 
H and pancreatic /3-cells [§,0]. The understanding of mechanisms underlying bursting is 
based on the fundamental work of Rinzel [II], accordingly to whom the neuron dynamics 
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implies two different time scales and it can be generally represented in the following form 

X = F(X,Y) 

(1) 

Y = G(X, Y) 

Here, the vector X describes the fast variables generating spikes and Y the slow ones con- 
tributing a modulation. The fast variables are associated with the membrane potential, 
Na + and K + conductance and other variables evolving on time scales comparable with the 
duration of the spike, whereas the slow variables refer to the concentration and conductance 
of Ca + ions Hl2| . Bursting arises as the evolution of the slow variables switches the fast 



dynamics between steady state and oscillatory dynamics. Thus a feedback process from Y 
confers to X the feature of a bursting pattern consisting of a fast spiking regime riding on 
a slow modulation [|12|,[n|. On the other hand, spiking phenomena are observed and easily 
modeled in other physical areas. Spike emission occurs in multi-mode class B lasers (solid 
state and semiconductor) operating close to threshold. This type of emission does not occur 



in a single- mode class B laser as CO2 , but it can be induced by a suitable feedback [ H| . 

Here we present evidence of the autonomous passage from a continuos train of irregular 
spikes to periodic bursts of regularly spaced spikes in a laser system. Precisely, we take as the 



fast system a laser undergoing homoclinic chaos of the Shilnikov type [|14|| . In this system, 
the presence of non-autonomous bursting induced by an external modulation has, already 
been reported |1| , with the irregular timing of the homoclinic spikes being regularized into 
a synchronized sequence of pulses. We consider such a behavior as generic for all spiking 
systems based on an activator- inhibitor competition (as e.g. biological clocks). Indeed the 
close approach to a saddle point, peculiar of homoclinic chaos, provides a local slowing down 
of the few variables (3-dimensional subset of X) describing the inflow to, and the outflow 
from, that unstable region. This means that all the other variables (complementary subset 
of X) do not play a crucial role around the singularity, whence such a generic homoclinic 
behavior applies to many different models which have in common the presence of a saddle 
point. The core features around the saddle point require three coupled variables; however, 
in order to recover the global aspects of phase-space orbits away from the saddle point, we 



will use a 6-dimensional model as reported later in Eqs.(2). Motivated by the fact that 
bursting in biological systems does not require an external forcing, we explore the feasibility 
of autonomous bursting in the laser. This is achieved by a secondary feedback loop, where 
a fraction of the output is fed back to a loss modulator after passage through a low pass 
filter. The slow dynamics thus provides the mechanism for bursting the laser intensity. 

The experiment has been performed on a single mode CO2 laser with an intracavity loss 
modulator; the detected output intensity is fed back in order to control the cavity losses 
(Fig. 1). The DC level and amplitude of the feedback signal are controlled by adjusting the 
bias (Bo) and the gain (R) of the amplifier. We set these two parameters so that the laser 
intensity displays a regime of Shilnikov chaos (the saddle point being a saddle focus [0). 
The laser intensity, after a time interval where it practically approaches the zero intensity, 
displays a large spike followed by a damped train of fast oscillations and a successive train 
of growing oscillations (Fig. 2). Damped and growing trains represent respectively the 
approach to, and the escape from, the saddle focus S from where the trajectory rapidly 
returns to zero and then starts a new orbit. In the chaotic region, the escape time from 
the saddle focus has different duration from pulse to pulse. The power spectrum shows a 
broadened peak at the average repetition frequency of the pulses, v n ~ 2 kHz. The non 
vanishing frequency content below v n suggests that it is possible to enhance one of those 
frequencies by feeding back a filtered fraction of the output. Thus, we filter the detected 
signal by a low pass filter with a cut off frequency (is c ) lower than v n and then re-inject it 
through a secondary feedback loop. As a result, taking also into account the phase response 
of the filter, we select a bursting frequency somewhat lower than, but proportional to, v c . 
The other low [y < u c ) frequency components are efficiently suppressed. This way, we 
obtain clusters of laser pulses periodically spaced at the bursting frequency. We have used 
"Rockland mod. 452" filters with a slope of 40 dB/decade. It was experimentally tested 
that a single filter is not sufficient; we were obliged to cascade two filters. In fact, we will 
see later that the numerical model requires at least three poles. The secondary feedback 
loop contains also a variable gain amplifier, AC coupled in order not to alter the primary 
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feedback loop setting. The filter yields bursting sequences with adjustable periods, reset by 
changing v c (Fig. 3). 

The reported experimental evidence on a single laser system should not mask the coop- 
erative character of the phenomenon here investigated. This character is demonstrated by 
the fact that such a feedback is dynamically equivalent to the coupling of two independent 
systems. For this purpose, we have measured the correlation time of the chaotic signal, 
which results to be around 100 fj,s, consistently lower than the average interspike interval 
[y~ x ~ 500/is). Once we select a feedback frequency around 500 Hz, well below u n , we are 
assured that the feedback onset time (~ 2ms) is much longer than the decorrelation time. 
Hence, our feedback is equivalent to the cooperative interaction of independent lasers. 

The dynamics is modeled by coupling the laser dynamics, described by a six- dimensional 
state vector X {xi : i = 1 — 6} JT7J, to the slow dynamics of a three-dimensional filter Y 
{.'/, : 3 = 1 " 3} : 

x\ = k xi(x 2 — 1 — k\ sin 2 (x 6 )) 

x 2 = — 71^2 - 2k x 2 x 1 + gx 3 + x 4 +p 

x 3 = -1ix 3 + x 5 + gx 2 + p 

Xi = -72^4 + 9x5 + zx 2 + zpo 

x 5 = -72X5 + zx 3 + gxi + zp (2) 
x 6 = -(3(x 6 + B - R(xi - (f)(y 3 - y 3 ))) 

m = -0(2/1 - xi) 

2/2 = -0 (2/2 - 2/i) 
2/3 = -0 (2/3 - 2/2) • 

In these equations the variable x\ represents the normalized laser intensity, x 2 and x 3 are 
respectively proportional to the population difference and sum of the two molecular levels 
resonant with the field, £4 and £5 represent the difference and sum of the populations of 
the rotational manifolds, which supply excitation energy to the resonant transition. The 
variable Xq is proportional to the feedback voltage and it affects the cavity loss parameter 
via the term &isin 2 (:r 6 ). The time is rescaled according to t = t ■ 7 • 10 5 . 71 ,72 , g and 
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are decay rates, po is the pump parameter, z is the number of sublevels contained in the 
rotational manifolds. In the Xq equation B and R are respectively the bias and gain of the 
primary feedback amplifier and <fi is the gain of the secondary feedback amplifier. The slow 
dynamics corresponds to a low pass filter with a third order pole at 9 = 2ni/ c , therefore it is 
modeled by the cascade of three linear integrators. The numerical values of the parameters 
have been chosen so that the free running behavior of the laser displays Shilnikov chaos 
(Table I). 

Fig. 4 shows the model behavior at the onset of the stable bursting corresponding to 
a perturbation of 2%. It is important to note that the use of a single pole filter does not 
provide bursting because the frequency component corresponding to the average repetition 
frequency of the chaotic pulses [y n ~ 2 kHz) is not sufficiently rejected. In this case, the 
action of the secondary feedback results in a modulation at the average repetition frequency. 

According to the general bursting model of Eq.(l), the filter variables represent the slow 
subsystem Y . In particular, the Y equations are similar to those used in the Hindmarsh and 
Rose model for neuronal bursting [T8| . 



In our case the bursting phenomenon is related to homoclinic chaos, which requires a fine 
parameter tuning in order to set the orbits around the saddle focus P^B}4T7[. As the control 
parameter Bo is moved above or below the homoclinic value, the dynamical system enters 
respectively a periodic behavior or it goes to a meta-stable fixed point (excitable system). 
This is clearly shown by a ± 1% stepwise change of B as shown in Fig. 5-a,b. Furthemore, 
in the case of a positive perturbation, the repetition frequency v T of the periodic spike train 
is a monotonic function of the perturbation amplitude, up to 5% (Fig. 5-c). 

Such an evidence demonstrates a fact that may be relevant for neurodynamics, that is, 
noise is not necessary to induce a bursting phenomenon. We have checked the independence 
of autonomous bursting from noise by applying to the laser pump an additive white noise up 
to 10 % of the pump's nominal value, without observing modifications. This demonstrates 
that autonomous bursting does not belong to the class of noise induced phenomena as e.g. 
stochastic resonance |1| or coherence resonance []20 . 



To make comparison with neurophysiological behavior, some real neurons do not show 
bursting when the synaptic connections with other cells are blocked. In such a condition, a 
single neuron typically shows chaotic spiking activity. Examples are provided by pancreatic 
/3-cells [^T| or by neurons of the pyloric Central Pattern Generator in the stomatogastric 
ganglion of crustaceans pifl . In these cases the appearance of bursting is related to the 
synaptic interaction among the cells [p^ -|25|]- This behavior can be interpreted in the frame- 
work of our model if we take the secondary feedback term as inducing a cooperative coupling 
among identical coupled lasers. As for the frequency v r of the regular spiking within each 
burst, Fig. 5-c shows that it depends sensitively upon the Bo perturbation. This suggests 
that an accurate control of Calcium concentration might be the neurodynamic tool for a 
reliable and sensitive coding of different stimuli, without need to recur to noise [p|, [26|l . 

In conclusion, the chaotic spiking regime peculiar of homoclinic chaos has been converted 
to an autonomous bursting regime by means of a feedback loop acting on a long temporal 
scale as compared to the average repetition time of the chaotic pulses. The observed regimes 
display characteristics similar to those observed in neurons and are reproduced by 

a model where the fast dynamics is modulated by a slow one through a secondary feedback 
loop; the frequency of the modulation sets the bursting period, while its amplitude sets the 
spike repetition frequency within each burst. 

We acknowledge partial support from the European Contract No. HPRN-CT-2000-00158. 
A. D. G. is supported by the European Contract No. PSS 1043. 



6 



References 



[1] J. Keener and J. Sneyd, Mathematical Physiology, (Springer- Verlag, New York, 1998). 

[2] B. Hille, Ionic Channels of Excitable Membranes, (Sinauer Associate Inc. Publishers, Sunder- 
land, Massachusetts, 1992). 

[3] C. von der Malsburg, The correlation theory of brain function, Internal Report 81-2 (Departe- 
ment of Neurobiology, Max Planck Institute for Biophysical Chemestry, Goettingen, 1981); 

CM. Gray, P. Koenig, A.K. Engel and W. Singer, Nature (London) 338, 334 (1989); 

S. Grossberg and D. Somers, Neural Networks 4, 453-466 (1991). 

[4] M. J. Berridge, P. E. Rapp, J. Exp. Biol. 81, 217 (1979). 

[5] P. E. Rapp, M. J. Berridge, J. Theoret. Biol. 66, 497 (1977). 

[6] J. L. Hudson, M. Hart, D. Marinko, J. Chem. Phys. 71, 1601 (1979). 

[7] M. Deschenes, J. P. Roy, M. Steriade, Brain Res. 239, 289 (1982); 
M. Steriade, Neuroscience 102, 243 (2000). 

[8] R. K. S. Wong, D. A. Prince, J. Neurophysiology 45, 86, (1981); 

F.K. Skinner, L. Zhang, J.L. Perez Velasquez and P.L. Carlen, J. Neurophysiology 81, 1274 
(1999). 

[9] I.C. Atwater, C.M. Dawson, A. Scott, G. Eddlestone and E. Rojas, J. of Horm. Metabol. Res. 
10 (suppl.), 100 (1980). 

[10] T.R. Chay and J. Keizer, Biophysical Journal. 42, 181 (1983). 

[11] J. Rinzel in Mathematical topics in population biology, Morphogenesis, and N euro sciences, 
Lecture Notes in Biomathematics, Vol. 71, edited by E. Teramoto and M. Yamaguti, (Springer- 



7 



Verlag, Berlin, 1987). 

[12] X.J. Wang, J. Rinzel in Brain Theory and Neural Network, edited by M.A. Arbib (The MIT 
press, Cambridge MA, 1995). 

[13] E.M. Izhikevich, Int. J. of Bifurcation and Chaos, 10, 1171 (2000). 

[14] F.T. Arecchi, R. Meucci W. Gadomski, Phys. Rev. Lett. 58, 2205(1987). 

[15] E. Allaria, F.T. Arecchi, A. Di Garbo and R. Meucci, Phys. Rev. Lett., 861, 791 (2001). 

[16] L.P Sil'nikov, Sov.Math.Dokl. 9, 163 (1965). 

[17] A.N. Pisarchik, R. Meucci and F.T. Arecchi, Phys Rev. E62, 8823 (2000); 
A.N. Pisarchik, R. Meucci and F.T. Arecchi, Eur. Phys. J. D13, 385 (2001). 

[18] J.L. Hindmarsh and R.M. Rose, Proc.R. Soc. London, Ser. B 221, 87 (1984). 

[19] L. Gammaitoni, P. Haenggi, P. Jung and F. Marchesoni, Rev. Mod. Phys. 70, 223 (1998). 

[20] A.S. Pikovsky and J. Kurths, Phys. Rev. Lett. 78, 775 (1997). 

[21] A. Sherman, J. Rinzel and J. Keizer, Biophysical J. 54, 411 (1988). 

[22] M.I. Rabinovich et al. IEEE Trans. Circuits Syst.-I 44, 997 (1997). 

[23] R.C. Elson, A.I. Selverston, R. Huerta, N.F. Rulkov, M.I. Rabinovich and H.D.I. Abarbanel, 
Phys. Rev. Lett. 81, 5692 (1998). 

[24] A. Sziics, P. Varona, A. R. Valkovskii, H.D.I. Abarbanel, M.I. Rabinovich and A.I. Selverston, 
NeuroReport 11, 563 (2000). 

[25] N.F. Rulkov, Phys. Rev. Lett. 86, 183 (2001). 

[26] V. Parpura and P.G. Haydon, Proc. Natl. Acad. Sci. USA 97, 8629 (2000). 



8 



Figure captions 
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FIG. 1. Experimental setup consisting of a CO2 laser with a primary feedback loop imposing 
a regime of homoclinic chaos and a secondary slow feedback loop (dashed block) which produces 
bursting. LT, laser tube; EOM, electro-optic modulator; D, HgCdTe detector; LPF, low pass filter; 
VGA variable gain amplifier (AC coupled); R and Bo respectively gain and bias of the amplifier in 
the primary feedback loop. 
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FIG. 2. a) Experimental time series of the laser output intensity in the homoclinic chaos regime 
(zero intensity level indicated), b) Zoom of a single pulse showing the approximate location of 
the saddle focus S at the start of the growing oscillation; notice that the preceding large damped 
oscillation occurs away from S and it is not part of the local dynamics at S. 
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FIG. 3. Experimental time series of bursting obtained with the Rockland filter for three different 
values of the cut off frequency a) v c = 100 Hz, h)u c = 300 Hz, c) v c = 600 Hz. The low frequency 
sinusoids are trains of the filtered feedback signal reported on a different scale; in fact the feedback 
amplitude is 2 % of the spike amplitude. The effective bursting frequencies are respectively a) 80 
Hz, b) 200 Hz and c) 350 Hz, that is, proportional to, even though lower than, the corresponding 
cutoff v c . 
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FIG. 4. Numerical results corresponding to the three cases reported in Fig. 3. a) 9 = 4.5 • 10 
(y c = 100 Hz), b) 9 = 1.3 • 10~ 3 {u c = 300 Hz), c) 9 = 2.7 • 10" 3 {v c = 600 Hz). 
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FIG. 5. Stepwise increase a) and decrease b) of control parameter Bo by ± 1% (dashed line) 
brings the system from homoclinic behavior to respectively periodic or excitable behavior, c) In the 
case a) the frequency repetition v r of the spikes within the regular train increases monotonically 
up to a ABq/Bq ~ 5 % ; above this value, there is a saturation. 
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Tablel 

Parameter values used in the simulations. 
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